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ABSTRACT 

We present integral field spectroscopy of the nebular line emission in a sample of 9 
brightest cluster galaxies (BCGs). The sample was chosen to probe both cooling flow 
and non-cooling flow clusters, as well as a range of cluster X-ray luminosities. The line 
emission morphology and velocity gradients suggest a great diversity in the properties 
of the line emitting gas. While some BGCs show evidence for filamentary or patchy 
emission (Abell 1060, Abell 1668 and MKW3s), others have extended emission (Abell 
1204, Abell 2199), while still others have centrally concentrated emission (Abell 2052). 
We examine diagnostic line ratios to determine the dominant ionization mechanisms in 
each galaxy. Most of the galaxies show regions with AGN-like spectra, however for two 
BCGs, Abell 1060 and Abell 1204, the emission line diagnostics suggest regions which 
can be described by the emission from young stellar populations. The diversity of 
emission line properties in our sample of BCGs suggests that the emission mechanism 
is not universal, with different ionization processes dominating different systems. Given 
this diversity, there is no evidence for a clear distinction of the emission line properties 
between cooling flow and non-cooling flow BCGs. It is not always cooling flow BCGs 
which show emission (or young stellar populations), and non-cooling flow BCGs which 
do not. 

Key virords: galaxies: clusters, cooling flows:general - galaxies: evolution - stars: 
formation - galaxies: stellar content 



1 INTRODUCTION 

The Brightest Cluster Galaxy (BOG) is a giant elliptical, 
often a cD, and is typically located at the center of the 
cluster's gravitational potential. The formation and evolu- 
tionary history of these large galaxies has remained an area 
of active research. Cl ose pairs and mult iple nuclei are com- 
mon in cD galaxies (|Laine et al.l |2003| ) , and hence it has 
been hypothesized that me rgers and galactic cannibalism 
jOstriker fc Hausman|[l977l ) drive their formation. 

However, cDs are often found at the peak of the cluster's 
X-ray surface brightness, where the hot intracluster gas can 
potentially cool to form molecular clouds, therefore an addi- 
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tional stellar or gaseous component, arising via the cooling 
intracluster medium, may be expected . Often BCGs in cool- 
ing f l ow clusters show Hq emission (I von der Linden et al.l 
I2OO7I : ICrawford et al.llT999l ). an exciting find as it is a sig- 
nature of current or recent activity in their cores. Diverse 
and dramatic morpho logies of the Hq em ission, such as long 
tails of emitting gas (iFabia n et al.ll200ll). a s well as highly 
concentrated emission ( Donahue et al and even more 

filam entary structures ( Conselice et al.lT2002l : iBlanton et al.l 
l200lD have been observed. 

lEdwards etahl (|2007l ) found that it is in the sample of 
cooling fiow BCGs, compared to the sample of BCGs as a 
whole, where there exists an increased likelihood of Ha emis- 
sion (found in ~70% of BCGs in cooling flow clusters, but 
only ~15% of BCGs in the sample as a whole). This implies 
that the cooling fiow status of the cluster is an important fac- 
tor for activity at the center. Furthermore, it suggests that 
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within cooling flow BCGs, cool molecular clouds, warm ion- 
ized hydrogen, and the cooling intracluster me dium (ICM) 
are related. This is consistent with the finding o f lEdee et alJ 
who showed that whenever there are detections of 
molecular hydrogen in a co oling flow cluster, ther e are also 
detections of Ha emission. ICrawford et al.l (|l999l ) obtained 
optical spectra of 256 X-ray selected BCGs and found that 
27% have emission lines and associate most of these with 
cooling flow clusters. 

In cases where the structure of the ionized gas is well 
studied, complex and irregular morphologies are often spa- 
tially correl ated with emiss i on fea tures in other wavelength 
regimes. As lPonahue et al.l ()200Cl ) found for Abell 2597 and 
PKS-0745-19, the ionized gas is filamentary, and similar in 
extent and structure to the molecular gas emitting in H2 
1-0 S(l) line. Other examples of well known Ha s tructures 
include Abell 1795, which shows a huge tail ( Cowie et al.l 
1 19831 ). as well as NGC 4696 in Cent aurus (Crawford et al . 
2005b) and NGC 1275 in Perseus (|Conselice et all I2OO2I : 
iFabian et al.|[200^ . Each of the last three show structures 
spanning tens of kiloparsecs and harboring the same features 
as seen for the X-ray emitting gas (Fabian et al. 2001; Craw- 
ford et al. 2005a,b). Furthermore, the CO (2-1) and Ha emis- 
sions are found to be associa ted in NGC 1275 (jSalome et al.l 
l2006h . lDonahue et all (|2000l ') found the H2 1-0 S(l) line emis- 
sion peak of NGC 1275 to be co-spatial with the peak of Ha 
emission, although confined to the very center of the BCG. 

Observations from the Chandra telescope have shown 
regions of star formation that are associated with bright 
lumps and filaments of gas whose radiative cooling times are 
short. In Abell 1795 for example, an excess of blue light is 
detected from the underlying central dom inant galaxy, sug- 
gesting a population of hot young stars (iMcNamara et al.l 
Il996l : lOegerle et ahlbOOll : iBildfell et al.ll2008l ). Additionallv. 
many rotationally excited transitions of CO been detected, 
as well as the O VI line, which indicates that gases at a 
range of temperatures exist (20 K and 100 000 K, respec- 
tively). Nevertheless, the amount of H2 emission is incon- 
sistent with gas cooling directly from X-ray temperatures 
into cool clouds. Cooling fiow models predict more cooled 
gas than is observed (Bohringer et al. 2001; Peterson et al. 
2003). Possibly, the mass is deposited into molecular clouds 
which are then reheated by one of s everal processes - ho t 
stars, shocks, or AGN, for example jWilman et al.|[2002h . 
and only a small fraction of the cooled gas is detected. The 
ICM and radio sources often appear to be interacting (Cro- 
ton et al. 2006) and these interactions may form cavities, or 
bubbles, in the surface brightness of the X-ray gas which 
move buoyantly through the ICM. In some systems, the 
bubbles carry enough energy to be able to balance the ra- 
diative losses e merging from the center of the clusters in 
the X-ray band ([P izzolato fc Sokei| l2005l : iBirzan et al]|2004 
iMcNamara fc Nul scn 2007). These AGNs are also potential 
sinks for the cooling gas and can contribute to the ionization 
of hydrogen gas clouds. 

In an attempt to distinguish between these scenarios, 
this paper examines detailed maps of the line intensity and 
morphology of the central regions of 9 BCGs. The line emis- 
sion properties are obtained using the integral field spectr- 
graphs. 

The objective is to identify and characterize the current 
and recent activity in BCGs. In order to develop non-biased 



conclusions, it is important to study the high X-ray luminos- 
ity systems, those of lower luminosity, and those in cooling 
fiows and non-cooling flows. The star formation rate (SFR) 
is calculated for young stellar populations and compared 
to the cooling flow mass deposition rate, which is calculated 
from the X-ray luminosity within the cooling radius (|Fabianl 
Il994h . Careful attention is paid to signatures associated with 
processes such as the cooling flow phenomenon, AGN activ- 
ity and galaxy-galaxy interactions. 

The paper is organized as follows. In section [51 we in- 
troduce our galaxy sample selection and outline our obser- 
vations. In section |3] we report our results: the emission 
line morphology and velocity maps, followed by an inves- 
tigation into the emission mechanism. In section |4l we con- 
sider the impact of our results on various galaxy and galaxy 
group formation hypotheses. Throughout o ur discus si on we 
will compar e our resul ts to those of Wil man et al.l (I2OO6I ) 
and Hatch et al.l (|2007l ) who have recently presented inte- 
gral held spectroscopy of the central regions of cooling flow 
BCGs. Unless otherwise stated our analysis assumes the val- 
ues 57ni=0.3 for the matter density parameter, 57a=0.7 for 
the cosmological constant, and _ffo=70 kms~^ Mpc~^ for the 
Hubble parameter, hx refers to the bolometric X-ray lu- 
minosity throughout. We will often refer to the thesis of 
I Edwards! (|2007h . within which additional detailed informa- 
tion can be foun43- 



2 SAMPLE SELECTION AND DATA 
REDUCTION 

Our nine targets are chosen such that line emission in the 
center of the galaxy is plausible. All BCGs are within 50 kpc 
of the X-ray center and either previous reports of emission 
lines exist (most notably from the catalogue of Crawford et 
al. 1999) , or the cluster properties are often associated with 
emission lines, i.e. a cooling flow (C F) is centered on th e 
BCG and radio emission is detected l|Edwards et al.l [20071 ). 
The images and spectra are observed using the integral field 
spectrographs on the Gemini Telescopes (CMOS IFU) and 
the Wilham Herschel Telescope (OASIS on the WHT). 

The majority of BCGs which have previously been ob- 
served with integral field units are in cooling flow clusters 
(Conselice et al. 2001; Crawford et al. 2005b; Wilman et al. 
2006; Hatch et al. 2007). Here, we include non-cooling flow 
cluster BCGs. The sample is listed in Table [1] The first col- 
umn lists the name of the cluster, the second the name of the 
BCG, and the third column lists the cluster redshift (taken 
from the NASA Extragalactic Database, hereafter NED). 
The distance and angular scale, assuming ho=0.70, are listed 
in the fourth and fifth columns, respectively. The dimensions 
of the integral field unit FOV are given in column 6. The 
cooling flow status (obtained from the literature) is given 
in column 7 along with the mass deposition rate (MDR) in 
column 8. The values in these last two columns are accompa- 
nied by two important caveats. First, the cooling flow status 
is not uniformly determined throughout the literature. Of- 
ten a short cooling time will be used to define a cooling fiow, 

^ An electronic version of the thesis can found at Scientific Com- 
mons http:/ / en.scientificcommons. org/301 78658\ 
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Table 1. Cluster Properties 



Cluster 
Name 


BCG 
Name 




D 
Mpc 


kpc/" 


IFU FOV 
kpc X kpc 


CF 

Status 


MDR 
Mgyr-l 


lO^'^W 


Abell 1060 


NGC 


3311 


0.0126 


53 


0.25 


1.3 


X 


1.8 


NqI 




0.47 


Abell 1204 






0.1706 


640 


2.65 


13.5 


X 


18.9 


Yes2 




6.77±1.42 


Abell 1668 


IC 4130 


0.0634 


256 


1.17 


5.9 


X 


8.2 


No3 




1.59±0.24 


Ophiuclius 






0.0280 


116 


0.55 


2.7 


X 


3.8 


No4,Yes* 




> 4 


MKW3s 


NGC 


5920 


0.0450 


184 


0.85 


4.3 


X 


6.0 


Yes5 




2.68±0.29 


Abell 1651 






0.0849 


337 


1.51 


7.6 


X 


10.6 


Yes* 




8.25 


Abell 2052 


UGC 


9799 


0.0345 


146 


0.68 


4.4 


X 


6.0 


Yes2 




2.52±0.20 


Abell 2199 


NGC 


2199 


0.0310 


125 


0.59 


4.4 


X 


6.0 


Yes6 


1211 


3.65±0.15 


Cygnus-A 






0.0561 


227 


1.04 


7.7 


X 


10.8 


Yes'^ 


~250 





The X-ray luminosity and errors (where availabl e) are fromlEbel inc^ et a7] \l 994 ) and lEbelina et al\ \l99A ). Refere nces for cooling flow 
status and MDR calculations: ^ Exosat data from \Havak awa et al! t200i); ^ Chandra and XMM-New ton data from] O 'Pea et al\ \2004) 
and^ASC A data from. Salome & C ombes (2oM) and ^hitt {200C); RGS XMM- Newton data from [Peterson et al\ \200A)~H ;handra 
data from \ Johnstone et al\ hOOj) : ROSAT data from lRevnolds & Fabiari tl99dl} and ^ Suzaku data from \Fuiitaetal. 



however, a high mass deposition rate, or a central temper- 
ature drop will also be used to classify the cluster. Overall 
these methods generally achieve the same results in terms 
of classifying the clusters, but can lead to discrepancies for 
the less powerful systems. Second, the mass deposition rates 
calculated based on the classical cooling flow paradigm are a 
factor of 10 to 100 times larger than those computed assum- 
ing most of the gas does not cool beyond ~1 keV. The values 
in the table quoted from ROSAT are the erroneous classic 
values and should be taken as much higher than the actual 
mass drop out rates. Note the very large errors on all MDR 
values. The X-ray luminosity is in column 9 (see the table 
notes for references). Each BCG and its surrounding local 
environment is shown in the acquisition images of Figure [1] 
The spatial extent covered is drawn on the image and spans 
5" X 7" for the CMOS IFU, and 10.3" x 7.4" for OASIS. 

The CMOS IFU observations were completed between 
February and June, 2006. This was done in queue mode us- 
ing the two-slit conflguration to allow for the largest fleld of 
view. The filter and grating pairings were chosen so to ob- 
serve Hof {\rest6563 A) at the redshift of the cluster. Three 
BCGs, and their standard stars, were observed using OASIS 
on the nights of June 28-29, and July 2, 2005. The 22mm 
enlarger was used. The MR661 configuration was used to 
obtain Ha and the MR516 configuration was used to ob- 
tain II/3 at Arest4861 A. Each target was observed close to 
the zenith in order to reduce the effects of atmospheric ab- 
sorption. Table[2]gives the instrument configuration, grating 
central wavelength, rest wavelength coverage and integration 
time for each galaxy. 



pies at both sides of the image. These artifacts, were caused 
by imperfect sky subtraction of the fiat fleld frames dur- 
ing the standard reduction procedures. To make a cleaner 
continuum image for Abell 1060, the reduction process was 
slightly modifled. Instead of using the pipeline for this case 
we manually constructed an image of the fiat fleld using only 
the same wavelength range as the continuum coverage. We 
then normalized this flat field image and divided it into our 
continuum image directly. Although this did not result in a 
perfectly clean frame, it did make for a noticeable improve- 
ment: the bright and dark fringes on the right side of the 
continuum image in Abell 1060 had a fiux difference which 
was 10% of the average continuum level, which decreased to 
a 2% difference in the reworked image. The variation in the 
high and low continuum level (discounting fringes) stayed 
constant at ~35% in each case. These artifacts disappear in 
the continuum subtracted line images, as both the contin- 
uum, and line -I- continuum images contain the fringes. 

Table [2j column 7, includes a value for the red dening 
expected from the Milky Way (ISchlegel et al.1 Il998''l. This 
galactic extinction was removed using the IRAF task dered- 
den. Subsequently, dopcor and the known cluster redshifts 
(see Table [l| are applied to deredshift the spectra. 

Each hexagonal lenslet is about 0.2" and subsampled 
onto a rectangular grid of 0.1" pixels. Each pixel's spec- 
trum is then median averaged with the value of its 8 closest 
neighbours in order to increase the S/N and better match 
the seeing (typically, ~0.8" to 1"). The pixels are slightly 
subsampled with respect to the seeing, however our results 
rely on the analysis of even larger regions. 



2.1 CMOS IFU 

Several Gemini specific programs from the package gem- 
tools in the spectral analysis software IRAF were used 
to reduce and analyze the spectra in the standard fash- 
ion (bias subtraction, cosmic ray rejection, fiat fielding, 
wavelength calibration, sky subtraction, and atmospheric 
extinction-correction). Baseline standard stars were used to 
perform the absolute fiux calibration. The final datacubes of 
Abell 1651, and most notably Abell 1060, show bright rip- 



2.2 OASIS 

The XOasis reduction package was used for the bias subtrac- 
tion, spectral extraction, flat fielding, wavelength calibra- 
tion, as well as for a first cosmic ray subtraction and absolute 
fiux calibration. Further cosmic rays were removed by me- 
dian combination of the frames, except when there were less 
than three exposures taken. In this instance we used IRAF's 
imreplace to create a bad pixel mask and the crfix routine to 
interpolate across the cosmic ray ruined pixels. To account 



4 Edwards et al. 



/ 



Abell 1060 



* * 



Abell 1204 



^Ophiuchus" 



« 



MKW3s 



Abell 1668 
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Abell 2052 



Abell 2199 



Cygnus-A 



Figure 1. Image of the cluster centers and the IFU position. The CMOS IFU (5" X T' ) or OASIS (10.3" X 74") field of view is 
marked as a box centered on the BCG. Abell 1060 : Taken in the r-filter around the BCG, NGC 3311. The central dust patch is already 
noticeable in this image. Abell 1204: Taken in the i-filter. A trail of light exists in the direction from the BCG to the small galaxies just 
East of North. Abell 1668 : Taken in the r-filter around IG 4130, the BGG. Ophiuchus: Taken in the r-filter. A hint of extinction is seen 
in the northern comer of the IFU FOV. MKW3s : Taken in the r-filter around NGG 5920, the BGG. Abell 1651 : Taken in the i-filter. 
Abell 2052 : Taken in the DSS red filter around UGG 9799, the BCG. Abell 2199 : Taken in the Gemini i-filter around NGC 6166, the 
BGG. Cygnus-A: Taken in the HST F622W filter around the BCG. The images are scaled in arbitrary flux units. 



for accidental pollution created in the housing chamber due 
to an LED, the sky subtraction itself was preformed using 
IRAF after the flux calibration was completed in XOasis; 
the entire sky frame was subtracted from the entire object 
frame, as the signature of the LED was apparent in both 
frames (each with the same integration time). Abell 2199 
had data taken in the MR516 configuration spread over two 
nights. After sky subtraction, flux calibration and cosmic 
ray removal, the datasets were averaged together. 

The absolute atmospheric extinction, A(0) for the WHT 
during the summer months is 0.12, and we used the stan- 
dard value of E(B-V) of 3.1 for our calculations of the atmo- 
spheric extinction. We used the NAOMI AO corrector and 
the seeing was ~0.7", which is greater than the instrumen- 
tal spatial resolution (0.26"). Therefore, the surrounding 8 
pixels were added to increase the S/N and match the seeing. 
The corrections of galactic extinction and dereddening were 
done using the same methods as in the CMOS data. 



2.3 Line Measurements and Internal Extinction 

The spectral line characteristics of flux, central position, 
and line width (FWIfM), are mapped in Section 3. We 
fit multiple Gaussian profiles to the spectrum of each IFU 
pixel using the IRAF task deblend. For the nearby Ha and 
[N ll] AA 6548,6584 emission lines we find the best fit for all 
three lines simultaneously. The task works by taking user 
specified best guess parameters for the continuum level and 
location of the line centers. It then varies the parameters 
of the line model in order to improve the chi square, based 
on the Levenberg-Marquardt method. Errors quoted on fiux 
levels and central position are also taken from the deblend 
fit. They are based on the Poisson statistics of the model 
of the data and determined by fifty monte carlo simulations 
which are run to create random Gaussian noise. The dif- 
ference of the fitted flux to the model flux determines the 
error. 

For most galaxies we are able to measure emission 
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Table 2. Observational Data 



Observatory 


Cluster NarQC 


Configuration 


A 


ZA A 

A 


int 1 ime 

sec 


rj{D—\ }Gal 

mag 


tijiyti— V )int 

mag 


GS 


AbcU 1060 


R400+r 


6300 


5710-6830 


1640 


0.079 


0.18 


GN 


Abcll 1204 


R400+i 


7800 


5940-6980 


3720 


0.017 


03 


GN 


Abell 1668 


R400+r 


6300 


5560-6615 


3600 


0.032 


0.3 * 


GS 


Ophiuchus 


R400+r 


6300/6350 


5450-6850 


6000 


0.591 


0.3 * 


GN 


MKW3s 


R400+r 


6300/6350 


5450-6700 


5400 


0.035 


0.3 


GS 


Abell 1651 


R400+i 


7800/7850 


6460-7700 


6000 


0.027 


0.3 * 


WHT 


Abell 2052 


MR661 


6610 


6130-6740 


3600 


0.037 


0.22 5 


WHT 


Abell 2052 


MR516 


5160 


4700-5345 


1800 


0.037 


0.22 5 


WHT 


Abell 2199 


MR661 


6610 


6430-6770 


2400 


0.012 


0.10 5 


WHT 


Abell 2199 


MR516 


5160 


4700-5345 


4800 


0.012 


0.10 5 


WHT 


Cygnus-A 


MR661 


6610 


5915-6605 


3600 


0.381 


0.6 6 


WHT 


Cygnus-A 


MR516 


5160 


4695-5225 


3000 


0.381 


0.6 « 



° Taken from NED. 



References for internal extinction: \ Vasterbera et al. 1 tl99i) : ^Sadler & Gerhard tWSd): ^ There is no correction for internal extinction 
in this case since the Balm er decrement could n ot be calculated i n\Crcm}forT'et al\ \l99A ): * We use 0.3, average value for the BCGs in 
the \ Crawford et al\ 111 99^ ) sample; ^ Taken from [Crawford et al\ tl99!Ji) : and^ Average value from this data. 
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Figure 2. Extinction map of Cygnus-A. The image is a rep- 
resentation of E(B—V) derived from the Ha to H/3 line strengths, 
with the Ha emission overlain as contours. One pixel is r~^220pc 
across. 



lines in Ha, [N ii] AA 6548,6584 and [O iii] A 5007 
and the absorption line of NaD at 5890 A with pixel- 
to-pixel S/N>10. In some cases we could also detect 
[O l] AA 6300,6364, and [S ll] AA 6716,6734 in emission with 
S/N>5. For the GMOS individual spectra, the minimum 
detectable level, la is ~3x 10~^^ erg s"^cm~^ (For the OA- 
SIS spectra obtained with the MR661 configuration, la is 
~lxlO"^* ergs"^cm"^). For the OASIS spectra taken with 
the MR516 configuration, la is ~7x 10~^* erg s~^cm~^ and 
absorption lines of H/3, Fe l and Mgj can only reach a 5a de- 
tection by adding the spectra over several individual pixels. 
Unless otherwise stated, the continuum level in the red has 
been determined by taking a median of a ~100 A- wide region 
around 6880 A if observed using the R400+i configuration, 
and around 6450 A if observed in the R400+r configuration 
(exact continuum regions can be found in Edwards 2007). 
The rms for each spectrum was calculated in the same spec- 
tral windows used to make the continuum image. Average 
values for the signal-to-noise (S/N) ratio of the continuum 
near Ha are between 10 and 20 (between 8 and 12 for OA- 
SIS). For both of the GMOS IFU configurations, R400-fr 



and R400-|-i, the observed resolution is 3.3 A, an average of 
the measured FWHM of observed sky lines at ~6300 A. The 
configurations used for OASIS each give an observed FWHM 
of ~2.5 A for the standard stars. 

In order to calculate qualities such as metallicity and 
age, and to compare the emission between galaxies, a 
knowledge of the internal reddening of each galaxy is im- 
portant. Some elli ptical galaxies are known to be dusty 
ISadler &: Gerhardlll985) and this dust can be uniform, fil- 



amentary, or patchy ( Laine et al.ll200^ . Our data allow a 
map of the internal extinction only for Cygnus-A and we 
show the work for deriving the pixel to pixel extinction for 
this case in the following paragraph. For all other galaxies, 
we only have access to an integrated value yielded from slit 
measurements obtained from the literature. This is by no 
means ideal, but as only the central few kpc are covered in 
the FOV, we resign to using this method. These values are 
listed in column 8 of Table [2] In some cases, no known val- 
ues of the internal extinction for a particular galaxy can be 
found and so we adopt a value of E(B— V)int ~ 0.3. This 
value is the average value for X-ra y selected BCGs w ith 
strong Ha emission lines found in iCrawford et all l|l999l ). 

For Cygnus-A, we align the images to the peak in the 
[O III] A 5007 line image to that in the [N ii] A 6584 line 
image. These two lines are chosen since they are both strong, 
and both are high ionization lines. The method is validated 
as the result is a good match to the continuum peak of both 
images. The extinction map is subsequently created using 
the following equation on a pixel by pixel basis: 



E(B-V),„t 



2.177 
-0.37 R 



(1) 



Here, R = 3.1, loHa/loU/i is 2.85, the theoretical ra- 
tio for Case B recombination, and lua and Ih/3 are the 
observed values (Osterbrock & Ferland 2006; Kauffmann 
et al. 2003). The theoretical IoHa/IoH/3 ratio of 2.85 may 
not be the ideal value to use for known Seyfert galaxies, 
like Cygnus-A, but the actual value is debated. It is of- 
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Table 3. Summary of IFU Emission Line Properties 



Cluster 


Lines 


SFR 


Ageoid 


Mass poid 


Name 














Mgyr-i 


Gyr 


M© kpc-2 


A1060 


SF 


0.02±0.001 


10±6 


2±1 


AI204I 


SF & AGN 


7.0±0.4 






A1668 


AGN 




7±3 


3±2 


OphiscG 


No2 




10±5 


150±80 


MKW3s 


AGN3 




>12 


>8460 


A1651 


No 








A2052 


AGN 




>4 


>30 


A2199 


AGN 




> 1 


> 3 



Cyg-A AGN 



Note 1: Possible companions near line emitting gas in A 1204- 
Note 2: There is a clump of gas, Object B, near to the BCG in 
Ophiuchus showing LINER or AGN-like emission line ratios. 
Note 3: Emission lines in MKWSs are blueshifted with respect to 
underlying cD population. 



ten assumed that the Ha emission in these systems is en- 
hanced due to coUisional processes, and several authors use a 
value of 3.1 (jGaskcU & Fcrl and 1984 ; Ostorbrock & Forland 
I2OO6I). although other values have also been determined 
(jSinette et al.lll99(]l . calculate a value of 3.4). The extinction 
map is presented in Figure [2] and shows very high values of 
E(B— V)int — 1.0 - 1.2 at the center where the Ha emission is 
maximal. Only pixel values with S/N>5 for both Ha and H/3 
are plotted. Although the values of extinction determined 
here may be slightly overestimated due to the choice of in- 
trinsic IoHq/Ioh/3 used, the average value in the NW emission 
peak from our map is 0.75± 0.2. This is not far from the inte- 
grated value of 0.69±0.04 (jOsterbrock fc F crland 2006) for 
Cygnus-A, also calculated using the Balmer decrement but 
using an intrinsic IoHa/IoH/3 ratio of 3.08. 



3 EMISSION LINE MORPHOLOGY, 
KINEMATICS AND DIAGNOSTICS 

We present the morphology, kinematic structure, and line 
diagnostics extracted from the continuum emission and line 
fluxes. Gaussian profiles are fit on a pixel-to-pixel basis, ex- 
cept for the BCGs in Abell 1668 and MKW3s, where single 
Gaussian fits are impossible and the flux is measured within 
a window around the position of the Ha line. The emis- 
sion line maps have all been constructed from continuum- 
subtracted measurements and no absorption correction is 
made. The region spectra shown are the integrated values 
of ~ 30 spectra and are labeled on the the Ha images. 

In terms of the morphology, we compare the results from 
the Ha and [N 11] A 6584 lines to those in the surrounding 
continuum for each system, and take note of any particular 
associated galaxy characteristics such as prominent dust fea- 
tures or nearby neighbou r s. Ma ps of [S 11] A 6716 -|- A 6734 
are available in lEdwardd (|2007l ). but not presented here as 
their morphology mirrors that seen in the [N 11] A 6584 maps. 
H/9, [O III] A 5007, and maps of the continuum near H/3 at 



Arest486lA (between 5025 and 5100A) are presented where 
available (Abell 2052 and Cygnus-A). Only pixels for which 
the flux measurement is >5a are used in the maps presented. 

To examine the kinematics, we show the velocity of the 
emission lines relative to the cluster radial velocity. When 
available, the NaD absorption line originating from the un- 
derlying galaxy is used for comparison. The velocity and 
FWHM maps also only show the pixels where the line flux 
measurement is >5a. 

To diagnose the origin of the emission, we present spec- 
tra for representative regions of the BCGs (usually at the 
location of the Ha or continuum peak). Generally, we do not 
have the H/3 and [O lll] lines re quired for constructing a BPT 
diagram (|Baldwin et al.lll981^ for each galaxy. However, we 
have plotted values of the [N 11] A 6584/Ha line ratio as a 
function of Ha luminosity for central regions in all the BCGs 
with emission lines. This flgure helps separate the ionization 
mechanism as a value of ([N 11] A 6584/Ha) > 0.63 is likely 
to be from an AGN or LINE R, rather than star formation 
l|Osterbrock fc Ferlandll2006l ). 

Before discussing the details of each system individually, 
we offer the reader a quick summary. 

• All cluster BCGs show smooth continuum emission, ex- 
cept in the case of Abell 1060. Although the continuum emis- 
sion is usually smooth, the morphology of the line emission 
is not uniform throughout the sample. Two BCGs show fil- 
amentary emission (in Abell 1668 and MKW3s). There are 
also two clusters for which the BCG harbours extended emis- 
sion (Abell 1204 and Abell 2199), as well as two in which the 
BCG line emission is condensed (Abell 2052 and Cygnus-A). 
In the cases of Abell 1060 (in the BCG) and Ophiuchus (just 
outside of the BCG) , there are large patches of dust and the 
line emission follows a similar morphology. Two BCGs show 
no evidence for line emitting gas (that in Abell 1651 and in 
Ophiuchus) . 

• The line kinematics usually vary smoothly in terms 
of velocity and line width. At times rotation is clearly 
present (Abell 1060 and Cygnus-A). Also, bulk motions 
(Abell 1204, Abell 1668, and Ophiuchus) and outflows 
(MKW3s) are observed. The subtraction of the absorption 
spectrum for MWK3s allows us to see a clear velocity shift 
of ~12 A (550kms~^) for the emission lines. 

• Except for the BCG in Ophiuchus and in Abell 1651, 
emission lines of Ha and [N 11] AA 6548,6584, are prominent 
(Figure [3|. For the three galaxies observed with OASIS, 
there also exist spectra around the H/3 line. The la noise 
level usually dominates any emission or absorption in H/3 
for Abell 2052 and Abell 2199. However, the [O ill] A 5007 
emission rises above the 5a level in central region (Region 2). 
Most of the galaxies show AGN-like spectra as summarized 
in Table [3] with very strong forbidden lines with respect 
to the Balmer emission. Figure [3] shows this for all galaxies 
with Ha emission. These lines also broaden in some regions 
of Abell 1204 and Cygnus-A. The BPT diagram of Cygnus-A 
confirms a Seyfert nucleus in this system. The lines in Ob- 
ject B of Ophiuchus are well described by ionization from 
a hard source, such as an AGN, but Object B is not the 
BCG. In the case of Abell 1060 and Abell 1204 for which 
the emission line ratios show very strong Ha emission with 
respect to the [N 11] A 6584 lines, we may be detecting a 
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Figure 3. Spectra of representative regions for each BCG. The spectra are the integrated values of individual spectra that 
make up the labeled regions identified in figures of Section 3. In the case of MKWSs the spectrum has been subtracted by the underlying 
galaxy spectrum, where absorption lines are shifted by -560±50 km s^^ . Blue spectra for the OASIS clusters are also shown. The size 
of the region covered for Abell 1060 is 175pc on either side, lllOpc for Abell 1204 (both regions), llOpc for Abell 1668, 330pc for 
Ophiuchus (BCG and Object B), 580pc for MKWSs, 650pc for Abell 1651, 710pc for Abell 2052, 630pc for Abell 2199, and 890pc for 
Cygnus-A. 
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Figure 4. Emission line mechanism. The [Nil] X6584/Ha ratio is shown. For NGC 3311 in Abell 1060 one pixel is ~675pc^, for 
the BCG of Abell 1204 one pixel is ~0.07kpc^, for IC 4130 in Abell 1668 one pixel is ^0.014 kpc^ , for Object B one pixel is ^3100 pc 2, 
for NGC 5920 in MKW3s individual pixels can not be plotted. The region size is r^7.29kpc ^. For UGC 9799 in Abell 2052 one pixel 
is r-^ 0.032 kpc^ , for NGC 6166 in Abell 2199 one pixel is ^0.023 kpc ^, and for Cygnus-A one pixel is ^0.0^8 kpc^ . A simplistic cut for 
AGN vs. SF is at y=0.63 (see also FiaurelM). 
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young stellar population. We determine the SFR for these 
last two. 

3.1 NGC 3311 in Abell 1060 

Morphology - Figure [7] shows the images of the Ha and 
[N ll] emission lines, as well as the continuum. Note that 
for this case the [OI] emission lines near 6300 A do not rise 
above the la noise level and therefore we use the median 
throughout a larger window (6300 and 6500 A) to build the 
continu um map. 

,Vasterberg et al.l (ligQll ) and iLaine et all l|2003h both 
note a large dust patch that corresponds in position to the 
obscuring feature going from the SW to the NE seen in the 
continuum image. In compa ring the HST I-ban d image of 
the center of the BCG from iLaine et al.l (|2003h . we notice 
the striking similarity both in extent and morphology be- 
tween the dust patch they identify, and the regions we iden- 
tify as strongly emitting in Ha. The Ha contours on the 
continuum image indicate that much of the Ha emission is 
confined within the dust, and therefore much Ha emission 
could be obscured. Overall, the bright regions in [N ll] fol- 
low those seen in the Ha emission image. The image of [S ll] 
A 6716 (not shown) follows the same overall morphology. 

Kinematics - Figure [H] shows maps of the relative ve- 
locity for the Ha and [N ll] A 6584 emitting gas. The ap- 
pearance of both emission lines is smooth, and we measure a 
velocity shear of 100±20kms~^ across ~1.0kpc. There is a 
clumpy distribution in the NaD velocity (clumps of -100 to 
-1-100 km s~^), originating from the underlying galaxy, but 
no shear. The emission line widths range from 130 kms~^ to 
200kms~^, close to the resolution, and have a clumpy dis- 
tribution (available in Edwards 2007). The velocity shear, 
low velocity values, and lack of structure in line widths are 
an indication of rotation of the line emitting gas. 

Emission Diagnostics - It is clear from the spec- 
trum presented in Figure [3] that none of the regions in 
NGC 3311 show ratios of significant AGN contamination. 
Figure m shows the line ratios of individual pixels and high- 
lights the area of Region 6, where the Ha emission is the 
strongest, and the line ratios are the lowest, falling well be- 
low the AGN/SF cut of 0.63. Based on the line ratios [N ii] A 
6584/Ha and also ([S ii] A 6716 -I- A 6731)/Ha, most regions 
confidently fall into the composite, or star formation section 
of the BPT diagram. The ([S ll] A 6716 + A 6731)/Ha are 
not shown here, but all regions have a ratio less than 0.25, 
much below 0.40, which is typical of an AGN or LINER. 
The image of the [N l l] A 6 584/Ha across the center of the 
galaxy, see lEdwardi ((2003), shows the central regions to 
have the lowest [N ll] A 6584/Ha ratios. Thus, it is fairly 
certain that the ionization is not a result of a hard radiation 
field supplied by an AGN or LINER. 

Star Formation Rate and Stellar Populations 
- Once the possibility of an AGN has been discarded, a 
common ionization mechanism explored is the effect of a 
young populatio n of hot stars, which can produce strong 
Balmer emission llCrawford et al. ll 19991 : iDonahue et al I2OOOI : 
IWilman et al.ll2006l ). Another possibility recently calculated 
by Ferland et al. 2008 is non-radiative heating by cosmic rays 
which can produce the observed molecul ar hydrogen lines (at 
least in the thin filaments of Perseus) . iFabian et al.l (|2008l ) 
have shown that star formation within the thin filaments of 
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Figure 5. Radial cut of the underlying population for 
NGC 3311 in Abell 1060. Based on the continuum flux near 
Ha. The x-axis shows the distance from the center of the IFU 
frame. The errorbars show the standard deviation of the con- 
tinuum values at the distance from the IFU center. There is no 
evidence for a significant trend with radius and most of the points 
are equal to each other. 



Perseus may be delayed by magnetic fields. However, at the 
larger distances of our clusters, it is unlikely that the ma- 
jority of the flux we observe is from such thin surrounding 
structures, and it is likely dominated by the bright emis- 
sion at the center of the BCG. Hence, for this paper, we 
explore the properties of a young population of stars that 
could excite the observed Ha emission. 

A lower li mit on the total SF R has been calculated using 
the method of iKennicuttI l| 19981 ). Excluding the dependence 
on the cosmology chosen, the errors in SFR are most sensi- 
tive to the error in Ha equivalent width measurement. The 
intensity of Ha emission can be diminished by the presence 
of Ha in absorption. The absorption comes not only from 
the same ionizing population of young stars, but also from 
the underlying and massive older stellar population of the 
galaxy. Although the absorption is below the la level in the 
individual pixel spectra, it can be removed in the integrated 
spectrum of the totality of Ha emitting pixels. We describe 
this presently. 

Absorption by the Underlying Population - The emis- 
sion appears only in the central 20"x35" (500pcx900pc for 
Abell 1060) of the IFU field of view, thus, a median of the un- 
derlying galaxy spectrum can be isolated by ignoring these 
central pixels with emission. We assume the older population 
has the same characteristics as (age and metallicity) , and is 
the dominant luminosity source for the overall luminosity 
profile of the optical image (that is, over the entire r filter). 
Therefore, the absorption should scale to this profile. Fig- 
ure[5]shows the average continuum flux in 5 pixel-wide radial 
bins where the errorbars show the standard deviation of con- 
tinuum values in the bin. The region 15-35 pixels (~400 to 
900 pc) from the center is plotted. The flux varies less than 
10%, from 3.1 to 3.5x10"^* erg cm"^ s"\A"\ showing no 
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(or only a weak) trend with radius. Therefore, we simply 
subtract the integrated absorption spectrum form the in- 
tegrated emission spectrum (scaled to the same projected 
physical area) . The average of both these spectra are shown 
in Figure [6] 

From Figure[6]we measure the integrated Ha absorption 
to be 3.2±0.4 x 10~^^ erg s~^ cm~^ per spectrum, on aver- 
age. With the absorption spectrum in hand, it is straightfor- 
ward to calculate an age and metallicity of the underlying 
stellar population by c omparing to the results of the popu- 
lation synthesis code of iMoUa fc Garcia- Varga3 (|2000l ). 

For Abell 1060 the NaD line is available. From the in- 
tegrated spectrum of non-Ha emitting pixels, an equivalent 
width of 4.5±0.4A, for the absorption is measured. Within 
our model, this strong absorption can only be fit with a 
super-solar metallicity {2Zq). An age of 1.0±0.6 x 10^" yr 
provides the best match. The error is based on the subset 
of models which fit the equivalent width within its mea- 
surement error. An estimate of the mass is made by scaling 
the continuum level (per solar mass) at 6400 A of the best 
model to the observed continuum level of the old popula- 
tion. For Abell 1060, 2.0±0.8 x 10* Mq kpc"^ are required 
in order to match the observations. Table |3] lists the age 
and mass of the older stellar population for this, and the 
other clusters where NaD, Fe I, Mgb, or H/3 absorption lines 
are available (we observe no absorption lines in Abell 1204, 
Abell 1651, and Cygnus A; details in the derivation of the 
old stella r popu lations for all the BCGs can be found in 
lEdwardd l\200ii )). Lower limits for the old population are 
quoted for Abell 1668 and Abell 2199 as low level Hq is 
observed in every pixel within the field of view. 

Star Formation Rate - We measure a total Ha fiux of 
5.0±0.3 X 10"^^ erg s"^ cm'^ or an Wo(Hq) = -106 A, 
in the absorption corrected integrated spectrum from 
all of the Ha emitting pixels. This results in a total 
SFR of 1.4±0.1 X lO-^Mgyr-^ The SFR density is 
1.2±0.1 X lO"^M0yr"Vc"^. The SFR derived is surely 
a lower limit as Figure [7] suggests, and a global measure 
of the dust absorption in the galaxy has already been ac- 
counted for. Assuming that most of the obscured emis- 
sion is behind this lane, a reasonable estimate on an up- 
per limit to the SFR would be to suppose that all of the 
obscured emission has the same intensity as in Region 6 
which is not behind the dust patch. Therefore, there should 
be no more than 1.7±0.1 x lO~^M0yr~^ in total, or, 
1.5±0.1 X 10^'^ M0 yr~Vc"^. This smaU amount of star 
formation found is less than t he one derived for cooling 
flow cluster BCGs. For example, iHicks fc Mushotzkvl (|2005l ) 
found rates of 0.2-219 Mr yr~^ derived from UV excess for 
cooling flow BCGs, and lEdwards et al] (|2007l ) found typical 
SFRs of 0.3 - 1.6M0yr"^ for emitting BCGs. The origi- 
nal MDR of 6M0yr"^ found for Abell 1060 would be be- 
tween 0.06 and 0.6 M0 yr~^ (1-2 orders of magnitude be- 
low the -RO^'^r value ). However, the Chandra observations 
(lYamasaki et al .I2OO2I ) show no evidence for any central tem- 
perature drop, and so no cooling flow is present in this sys- 
tem. Thus, the fact that we have measured only a very small 
amount of star formation activity we derive is consistent 
with this picture. 

Properties of the Young Population - We can further 
quantify the activity by estimating a metallicity, age and 
mass for the ionizing population. We start by calculating the 



m etallicity of the gas, usin g abundance ratios as described 
in lKewlev fc Dopital l|2002l ). and then we assume the metal- 
licity of the young stellar population is the same as for the 
gas. Although this may not be exact, it will help us to re- 
strict the solution for the age of the young stellar population 
using synthesis models. 

It is possible to make a rough estimate of the 
metallicity in the regions by assuming an average 
ionization parameter, and comparing the ratios of 
[N 11] A 6584/( [S II] A 6716 + X 6 7 31) an d [N 11] A 6584/Ha. 
Although, as Kewlev fc Dopital (I2OO2I ) strongly caution, 
both ratios depend strongly on the ionization parame- 
ter and this method is not precise. Assuming the ioniza- 
tion par ameter is between ave rage values of 5 x 10^ and 
2 X 10'' (|Kewlev fc Dopita|[200^ . the metallicities from the 
[N 11] A 6584/([S 11] A 6716 -I- A 6731) ratios are compared 
against those from [N 11] A 6584/Ha. The latter generally has 
two possible values, therefore constraints from the former en- 
able the correct point to be chosen, although large uncertain- 
ties exist for the [S 11] A 6731 line measurement. Considering 
the region within the central 0.60" x 0.70" (~175pc across), 
we derive a metallicity of 12-flog(0/H) — 9.3±0.1 (approx- 
imately t wice solar). This i s high , but the mass-metallicity 
relation of lTremonti et al] (|2004l ) shows higher metallicities 
in more massive galaxies. 

The age of the young stellar population is esti- 
mated using the evo lutionary synthesis code Starburst99 
(|Leitherer et all 1 19991 ). The Ha equivalent width is calcu- 
lated and matched to results from a run based on an instan- 
taneous burst of star formation using a Salpeter IMF with 
masses between 1 and 100 M0. We use a model based on 
super-solar metallicities (Z=2Z0) to best match the metal- 
licity (although a solar metallicity does not significantly af- 
fect the age). The best fit model, experimenting with both a 
mass cut of Mup = 100 M0 and 30 M0 suggests young ages 
(<10Myr) as expected for these ongoing bursts in which the 
youngest stellar populations dominate the spectral energy 
distribution. 

The Starburst99 code also lists the theoretical spectral 
energy distribution for a lxlO®M0 burst as a function of 
wavelength. We convert the model continuum luminosity 
to a flux and scale to the continuum level of the absorp- 
tion corrected emission spectrum, 6 x 10^^^ erg s~^ cm~^ 
(the 3cr errorbars are ±2 x 10~^^ erg s~^ cm~^). At lOMyr, 
this requires that the young population is 3±lxlO^M0. 
The masses are lower for younger ages and lower theoret- 
ical Mup, for example, at 5Myr with M„p=3O0 the mass is 
l.O±O.3xlO'^M0. 

These results show that if the Ha emission does indeed 
result from a young population of stars, it is consistent with 
a low level of star formation, which in turn is much less than 
that predicted by direct cooling from the ICM. 

3.2 The BCG in Abell 1204 

Morphology - Figure [9l shows the images of the contirmum 
subtracted Ha, [N 11] A 6584, and the continuum near the 
Ha emission line for the BCG. The line images share the 
position of the peak intensity with that of the continuum 
emission (including [S 11] A 6716 which is not shown). They 
also reveal the existence of a plume of bright Ha emission 
which extends from the central peak to the North, and then 
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Figure 6. Average spectra of the emitting and underlying population of NGC 3311 in Abell 1060. A median average 
of Ha emitting spectra are shown alongside the average spectrum for the underlying population in the regions surrounding the intense 
emission. The hashed regions are those affected by the chip gap and poorly subtracted sky lines. 
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Figure 7. NGC 3311 in Abell 1060. From left to right, the continuum subtracted Ha emission flux, the continuum subtracted 
[Nil] A 6584 emission flux, and the continuum near the Ha emission line. Both emission line maps follow a similar morphology overall 
and much of the emission is found within regions of strong extinction in the continuum. In the continuum image, a dust patch in the 
N-S direction is clearly visible. The vertical bands seen on the left and right sides of the figure are associated to fringes. The Ha emission 
is overlain as contours. Th 
One pixel is ^26 pc across. 
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Figure 8. Relative velocities for NGC 3311 in Abell 1060. Left Panel: Map of the Ha relative velocity for NGC 3311 in Abell 1060. 
Right Panel: Map of the [N ll] X6584 relative velocity for NGC 3311 in Abell 1060. Both maps are consistent with rotating gas. The Ha 
emission is overlain as contours and the scale is in units of kms~^. One pixel is ^26pc across. 
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Figure 9. BCG in Abell 1204. From left to right, the continuum subtracted Ha emission flux, the continuum subtracted [Nil] X 6584 
emission flux, and the continuum near the Ha emission line. The Eastern plume of emission extends towards several possible neighbour- 
ing galaxies seen on the acquisition image. The Ha emission is overlain as contours. The regions used in the analysis are represented 
as boxes. The images are in units of 10~^^ ergs~^ cm~^ . One pixel is ^270 pc across. 
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Figure 10. Maps of the kinematics for the BCG in Abell 1204. Left Panel: Map of the Ha relative velocity Right Panel: Map 
of the Ha FWHM. No obvious ordered motion is seen in the velocity map, however emission near the galaxy center is blueshifted with 
respect to the average velocity. The FWHM increases closer to regions of higher Ha to [N n] ratios and closer to the AGN. The Ha 
emission is overlain as contours and the scale is in units of km s~^ . One pixel is ^270 pc across. 



towards the East out to the edge of the image. This is in 
the direction towards several smaller galaxies seen on the 
acquisition image (Figure [TJ. There is what appears to be a 
smaller galaxy in the western corner of the continuum image 
which has no counter part in the emission line images. 

Kinematics - Figure [TOl shows maps of the relative ve- 
locity and FWHM for the Ha emitting gas (with the Ha 
emission overlain as contours). Regions of slight blueshift- 
ing, as well as redshifting are apparent. The physical scale 
is from — lOOkms"^ to +150kms~^, with the most nega- 
tive values echoing the structure of the Ha emission. There 
is no obvious ordered motion such as the rotation seen in 
NGC 3311 of Abell 1060. There is however some structure 
in the relative velocity map as the central emission is clearly 
blueshifted with respect to the rest of the emission (by up 
to 200kms~^). The mean FWHM across the field of view is 
130±50kms~^, but broadens to 570 ±70 km s"'^ at the cen- 
ter, in Region 4. The velocity differences are not as high 
seen in strong outflows, but this region hosts the highest 
[N ll] A 6584/Hq ratios, suggesting that the emission lines 
are ionized by an AGN (to be discussed further below). The 
line widths would be consistent with the central region being 
closer to the broad line region of the AGN or LINER. 

Emission Diagnostics - Comparing Figures [3] and |3] 
we see that [N ll] A 6584/Hq: luminosity ratio of star form- 
ing Region 1 is lower than line ratio of the AGN affected 



central Region 4. The hue ratios of [N ll] A 6584/Ha and 
also ([S ll] A 6716 + A 6731)/Hq suggest that most regions 
(Regions 1, 2, 7, 8, and 10) are likely a composite of emis- 
sion from young stellar populations as well as a LINER, 
while Regions 3, 4, and 5 are from LINER emission only. 
We cannot distinguish between a Seyfert and a LINER since 
the ratio of [O ill] A 5007/H/? was not me asured. Inte- 
grated along their slit, ICrawford et al.l (Il999l ) measured a 
ratio [O III] A 5007/H/3 > 1 and [N ll] A 6584/Hq = 1.4. 
These ratios will not be uniform throughout the central re- 
gions of the BCG. Nevertheless, their integrated value for 
[N ll] A 6584/Ha is between those of our regions. It is there- 
fore useful to keep in mind their result of [O ill] A 5007/H/3 
as an average value, possibly closer to a LINER. 

Star Formation Rate - It is not possible to calcu- 
late a SFR for Regions 3, 4, and 5 for example, as they 
are clearly dominated by AGN signatures. It is also pos- 
sible that the SFR calculated in the other regions are af- 
fected by the AGN, so in this way the rates presented 
should be interpreted as upper limits, contaminated by 
AGN ionization. No strong absorption lines are present, 
so there is no absorption correction to the integrated Ha 
equivalent width. A total SFR is derived from the com- 
bined spectrum of the brightest 255 Ha emitting pixels 
and yields a fiux value of 14±1 x 10~^^ erg s"'^ cm~^. This 
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Figure 11. IC 4130 in Abell 1668. From left to right, the continuum subtracted Ha emission flux, the continuum subtracted 
[N ll] A 6584 emission flux, and the continuum near the Ha emission line. In this case the emission has not been derived by Gaussian 
fits but instead includes the integrated flux between 6554 6572 A for Ha and between 6575 and 6593 A for [N nj. Both lines maps 
follow the same overall morphology, though the peak is displaced. The Ha emission is overlain as contours. The regions used in the 
analysis are represented as boxes. The images are in units of 10~^^ ergs~^ cm~^ A~^ . One pixel is ^120pc across. 
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Figure 12. Relative velocities for IC 4130 in Abell 1668. Left Panel: Map of the Ha relative velocity. Right Panel: Map of the 
NaD relative velocity. The line emission is not at rest with respect to the underlying galaxy. In addition, the gradient in velocity seen in 
the line emission is not apparent in the underlying galaxy .The scale is in units of kms~^ and is normalized to the redshift of IC 4130. 
Contours of the continuum emission near Ha are overlain and one pixel is ^120pc across. 



yields a SFR of 7.0± O.4M0yr"\ or a SFR density of 
2.3±0.1 X 10"* M© yr"^ pc'^ 

This upper limit is not inconsistent with the revised 
value of 50^30 Mq yr"^ for the MDR assuming the gas does 



not cool below O.lkeV found bv lO'Dea et al.r(|2008l ). as not 
all the molecular gas will convert to stars, and as the cur- 
rent estimates of MDRs are an order of magnitude helow the 
previously derived rates . Our total SFR is very close to that 
derived bv lO'Dea et al.l ( 2008|) from infrared luminosities in- 
side of a 12.2" aperture i Quillen et al.|[2008h . especially con- 
sidering our effective aperture is slightly smaller. Our results 
further suggest that although not the dominant source, part 
of the luminosity in the O'Dea measures may arise from the 
central AGN. 

We do not attempt to predict a valid age and mass es- 
timate for this source, as contamination from the AGN is 
important, and the SFR is only an upper limit. The inter- 
ested reader may refer to Edwards 2007 for an examination 
into the ages in select less contaminated regions. 



3.3 IC 4130 in Abell 1668 

Morphology - There are several pixels which are not well 
fit by single Gaussian profile fits for this target because of 
the lower signal, and here, the flux is simply added within 



a specific window. In this case. Figure [TT] shows the flux 
added between 6554 and 6572 A for Ha, and between 6575 
and 6593 A for [N ll] A 6584. The figure also shows the con- 
tinuum near Ha. Unlike the smooth elliptical distribution of 
the continuum image, both the Ha and [N ll] A 6584 flux 
images show much more patchy and filamentary structure 
throughout. The bright regions on the [N ll] A 6584 image 
correspond to bright regions on the Ha image, and the two 
emission lines share the same overall structure. However, the 
peak emission is displaced in the images. 

Kinematics - A map of the relative Ha line velocities 
is shown in Figure [121 It shows a clear gradient from positive 
values North of the center 100 kms~^), to negative values 
— lOOkms"^) at the South with the zeropoint near the 
center of the continuum. The low S/N values at the edge 
of the figure, as well as a few central patches where emis- 
sion lines could not be fit as single Gaussians are plotted 
as having zero relative velocity. Also shown in the figure is 
the velocity of the NaD lines relative to the average veloc- 
ity of IC 4130. This map traces the motion of the underly- 
ing galaxy, as opposed to the line emitting gas. It shows a 
smaller variation in magnitude (~ ±70kms~^), no smooth 
gradient, and all of the velocities are negative. This implies 
that the line emitting gas is not at rest with respect to the 
underlying galaxy. 
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Figure 13. BCG in Ophiuchus. From left to right, the continuum subtracted Ha emission flux with the continuum emission overplotted 
as contours, the continuum subtracted [Nil] X 6584 emission flux with the Ha emission overlain as contours, and the continuum near the 
Ha emission line also with the Ha emission overlain as contours. Both emission lines show similar morphology and both are displaced 
from the center of the BCG, at the location of an extinction patch. The regions used in the analysis are represented as boxes. Extinction 
is clearly seen in the Northern comer of the continuum image. The images are in units of 10~^^ ergs~^ cm~^ . One pixel is -^SOpc 
across. 



Emission Diagnostics - All regions show stronger 
[N ll] A6584 line fluxes than those for Hq as can be seen 
in Figures [3] and |4] Region 8 has the smallest ratio of all the 
regions, but the value is 1.7±0.2, including the measurement 
error from both of the lines. This is still well described by 
AGN or LINER emission. The facts that the [N ll] A6584/Ha 
ratio does vary across the image and with Ha luminosity, 
and that the scaled underlying spectrum has a lower con- 
tinuum than the emitting spectrum (although some amount 
of continuum is surely from the nuclear emission) suggest 
a young population may still exist. However, its lines are 
completely masked by the presence of the AGN or LINER, 
and therefore, it is not possible to use the Ha emission line 
to characterize a SFR or an age for the younger stellar pop- 
ulations. 

3.4 Ophiuchus 

Morphology - The continuum subtracted Ha flux image 
for the BCG is shown in Figure 1131 The continuum sub- 
tracted [N ll] A 6584 flux image is also shown, but the con- 
tours represent the Ha flux. There is no line emission at the 
center of the BCG, rather it is concentrated North of the 
continuum emission. This concentration of line emission is 
hereafter referred to as Object B. 

There is a drop in the continuum levels of the BCG to 
the North side, seen in Figure 1131 This is also apparent in 
Figure [T] where the dust is seen as a decrease in flux just 
to the side of the bright line-emitting Object B. Both these 
figures show a clear lack of emission at the same projected 
position as Object B, probably caused by dust. Since the 
dust feature appears to be very localized relative to the BCG 
and it is absorbing the galaxy continuum we conclude that 
the dust patch is in front of the BCG. 

The Nature of Object B - By comparing the average 
redshifted positions of the emission lines in Object B to the 
average velocity of the BCG NaD absorption line we measure 
a velocity difference of -|-600±50 kms~^ between Object B 
and the BCG, see Figure [3] This is well within the typical 
velocity dispersion of a massive galaxy cluster. There is also 
a difference of -1-150 km s~^ between the BCG and the cluster 
velocity (derived from the cluster redshift measurement of 




Figure 14. Map of the Ha relative velocity for Object B in 
Ophiuchus. The velocity scale here is with respect to the average 
velocity in Object B but this is offset from the BCG velocity by 
800 km s~^ suggesting Object B is in front of the BCG and falling 
onto it. The Ha flux is overlain as contours. The scale is in units 
ofkms~^ relative to the average Object B velocity. One pixel is 
r^60pc across. 



z=0.028 from lLahav et all (jl989l )l. though this difference is 
close to the error of the cluster radial velocity measurement. 

Three scenarios can explain the nature of Object B: 1) 
that it is infalling onto the BCG, 2) that it is a background 
galaxy, and 3) that it is being ejected from the BCG. The 
first scenario is supported by the absorption seen in the con- 
tinuum and acquisition images as well as the relative veloc- 
ity difference between Object B and the BCG. Scenario 2 is 
supported by the velocity difference but the dust absorption 
of the BCG in the region of Object B would have to be a 
coincidence. If Object B were being ejected away from us, 
the third scenario could produce the observed velocity dif- 
ference, but would not explain the dust absorption. Finally, 
if Object B were being ejected towards us, the dust feature 
would be explained, but the velocity difference would not. 

The hypothesis for which Object B is in front of the 
BCG and falling onto the BCG or spiraling around it most 
completely explains our observations. Figure [14] shows the 
Ha velocity map measured with respect to the average Ob- 
ject B velocity. This map could indicate some rotation of 
Object B (where the side near the BCG is rotating in our 
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Figure 15. NGC 5920 in MKW3s. From left to right, Ha + [N nj AA 6548,6584 o,nd the continuum near Ha. In this case the line 
strength is measured by adding all the flux within 6530 and 6590 A. An elongated emission feature that runs across the center of the 
galaxy can be observed. The Ha emission is overlain as contours. The regions used in the analysis are represented as boxes. The images 
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Figure 16. UGC 9799 in Abell 2052. From left to right, the continuum subtracted Ha emission flux, the continuum near the Ha 
emission line. Both maps show smooth, centrally condensed emission. The Ha emission is overlain as contours. The regions used in 
the analysis are represented as boxes. The images are in units of 10~^^ ergs~^ cm^^ A~^. One pixel is ^150pc across. 



direction) perhaps perturbed due to some tidal effect with 
the BCG. In this scenario, Object B could be a small galaxy 
(or a large molecular cloud) where the interaction with the 
BCG is compressing the gas resulting in part of the optical 
and X-ray emission observed. 

The integrated spectrum of Object B has 
[N ll] A6584/HQ! = 3.2±0.3, a value which is similar 
to those in its individual regions of (Figure |3| , and which 
indicates an AGN or LINER ionizing source. 

Furthermore, [O l] A6300/Hq = 0.3±0.05, and 
[S ll] A6716 + A6731/HQ = 0.9±0.3 are consistent with the 
AGN region of the BPT diagrams, though less clear. The 
([S ll] A6716 + A6731)/Ha for instance, varies in the differ- 
ent regions, from 0.72 to 1.26, which is nearer to but on the 
border of the AGN or LINER side of the BPT diagram. 



3.5 NGC 5920 in MKW3s 

Morphology - Figure [15] shows the image of the contin- 
uum subtracted Ha + [N ii] AA 6548,6584. Deep absorption 
lines from the underlying galaxy are shifted in velocity with 
respect to the emission lines. This superposition renders sin- 
gle Gaussian fits to the emission lines impossible. Therefore, 
the image shown is the addition of all flux between 6530 
and 6590 A. In contrast with smooth spherical continuum 
emission, the line emission shows a strong elongated feature 




Figure 17. Map of the Ha relative velocity for UGC 9799 
in Abell 2052. There exists a gradient in the velocity, however it 
is not smooth enough to be able to differentiate between rotation 
or an outflow. The Ha flux is overlain as contours. The scale is 
in units of kms~^, and values for pixels with bad Ha line fits 
are plotted as having zero relative velocity. One pixel is '^ISOpc 
across. 
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which crosses the center of the cD galaxy from the NE to 
the SW. 

Kinematics - In Figure [S] we present the spectrum of 
the central region (Region 2) for which, in this case, the ab- 
sorption spectrum is subtracted. Only here, do we subtract 
the absorption, as in this case as it is deep and shifted in 
velocity from the emission lines. We measure the absorption 
by building a spectrum of the underlying galaxy, as in Sec- 
tion 3.1. After the subtraction, the emission lines become 
easier to distinguish in the integrated spectrum of Region 2. 
They are clearly blueshifted from the underlying galaxy by 
— 560±50 kms^^. This is a known radio galaxy and the ve- 
locities are probably indicative of an outflow in our direction. 

Emission Diagnostics - Individual pixels are of much 
too low S/N to plot on Figure |4l we therefore plot the values 
of the (absorption corrected) integrated regions. The average 
ratios are also calculated for Regions 1, 3, 4, and 5, there 
are close to or larger than for Region 2. The values are 1.3, 
1.3, 2.6, and 1.7. All emitting regions are thus likely ionized 
by a hard source, AGN or LINER-like emission. 

3.6 The BCG in Abell 1651 

The integrated spectrum shown in Figure [3] has no emission 
lines but only Ha in absorption. Thus, although this cluster 
has a cooling core, there are neither signs of gas ionized by 
an AGN, nor by a population of hot young stars at the center 
of the BCG. 



3.7 UGC 9799 in Abell 2052 

Morphology - Figure shows the images of the contin- 
uum subtracted Ha flux, and the continuum near Hq (be- 
tween 6380 and 6430 A). These images show smoothly vary- 
ing emission and share the same peak location. We do not de- 
tect any H/3 emission above the la level of the noise, the up- 
per limit in the central Region 2 being 4x10"^^^ ergs~^cm~^. 
Lines of [N ii] A 6584, [S ii] A 6716 -I- A 6734 and 
[O III] A 5007, and the continuum around H/3 (which are 
not shown) display a similar morphology as Ha. 

Kinematics - The relative velocity map of Ha, Fig- 
ure [TT] has a range of — 250kms~^ to +150kms~^. The 
figure shows negative velocities to the South of the center 
of the emission, and more positive velocities North of the 
center of the line emitting region. The the line widths (map 
not shown) vary by 450-700 km s~^ across the galaxy, but 
show no clear structure. However, compared to cases with 
clear signs of rotation where the velocity gradient develops 
smoothly from positive to negative along an axis, in this 
case, it is less clear as to whether there is rotation or an 
outflow. 

Emission Diagnostics - Figures |3] and |4] show ratios 
characteristic of ionization due to AGN or LINER activity in 
most of the pixels. Regions 1, 2, and 6 have [O ill] A 5007 line 
fluxes of 8.5, 8.2, and 2.4xl0~^^ ergs~^cm~^, respectively. 
There are the only three regions where the line flux is greater 
than 5(T (2.0x 10~^^ ergs~^cm~^). Therefore, in these re- 
gions, the [O III] A 5007 emission is necessarily stronger than 
the H/3 emission and the [O lll] A 5007/H/3 ratio is ;> 6. This 
helps to constrain the harder emission source as a Seyfert 
nucleus, over a LINER ([O ill] A 5007/H/3 >3.2; also, see 



Crawford et al. 1999 who obtain an O ill] A 5007/H/3 ratio 
of 9.5). The ([S ll] A 6716 + [S ll] A 6731)/Ha line ra- 
tio is above 1.10 for all regions and well within the AGN 
side of the diagram. Therefore, it is not possible to use the 
Ha emission line to search for any young stellar population. 
However, observing the scatter in the individual pixels of the 
central Region 2 on Figure |4] does suggest that there may 
be 2 cases: 1) for Ha luminosities that increase as the ratio 
decreases, and hence are affected by star-forming regions, 
and 2) where the ratios are higher and the Ha luminosities 
are ~1. 

3.8 NGC 6166 in Abell 2199 

Morphology - Images of the continuum subtracted Ha flux 
along with the continuum subtracted [N ll] A 6584, and 
the continuum near Ha are shown in Figure 1181 The con- 
tinuum image which includes flux between 6650 and 6700 A, 
increases smoothly in brightness towards the center, without 
any prominent dust features. The peak of the Ha emission 
does coincide with the peak of the continuum flux, however, 
a second bright peak, not noticeable in any of the other line 
images is seen to the South. The Ha line flux continues to 
extend towards the East, a morphology that is mirrored in 
the [S ii] a 6716 + A 6734 emission line image (not shown), 
although in this case the [N ll] A 6584 line image differs in 
that it is much more condensed. The elliptical morphology 
of the continuum in the blue, between 5025 and 5100 A (not 
shown) is similar to that seen for the red continuum, though 
it is less condensed. The individual pixels of the blue config- 
uration show no strong H/3 or [O ill] A 5007 line emission. 

Kinematics - The relative velocity map of Ha emission 
is shown in Figure [T^ The velocities are calculated with re- 
spect to the rest frame of the central galaxy and plotted as 
such in the figure. However, one can see that clearly the emis- 
sion lines are not at rest with respect to the central galaxy 
(zflCG=0.03035±0.00003; de Vaucouleurs et al. 1991), and 
that the velocity difference within the line emitting gas itself 
goes from — 200kms~^ to -1-200 kms~^. The low level flux 
of the Eastern side of the image has a complex structure. It 
is the location of both the Region 6 with the largest veloc- 
ity, 700±50kms~^, and Region 5 with the lowest velocity, 
580±50kms~^. The pixels which have the most intense line 
emission are those with the lowest velocities, it could be that 
they are closer to us and hence subject to less extinction. 
If this were the case, one would then expect the H/3 line 
emission to be more smooth, however hardly any emission 
lines are seen in the image (only Region 8 has detectable H/3 
emission). However, the H/3 emission lines are quite possi- 
bly effected by absorption which cannot be seen at the S/N 
level of this data. In this case we would be seeing a curved jet 
with relative velocity ~500kms~^. A higher S/N H/3 image 
would help verify this scenario. 

Emission Diagnostics - The spectra in the red and 
blue are shown in Figure [S] The [N ll] A 6584/Ha ratios (see 
Figure |4)| for the different regions in this galaxy vary from 
1.7 to 10. The [O ill] A 5007/H/3 ratio, measurable only for 
Region 8, h as a value of 1.8 = fc0.5. This is consistent with the 
ratios from I Crawford et al.l |l99^ derived from a 6" slit. It 
is possible that stellar absorption affecting Ha and H/3 could 
be artificially raising the ratios, and hiding the effects of a 
young population. However, this would not be a significant 



IFU Observations of BCGs 17 





20 25 30 



5 10 15 20 25 30 




Figure 18. NGC 6166 in Abell 2199. From left to right, the continuum subtracted Ha emission flux, the continuum subtracted 
[N u] \ 6584 emission flux, and the continuum near the Hex emission line on the right. Neither emission line shows emission as 
condensed as seen in the continuum. The Ha emission is overlain as contours. The regions used in the analysis are represented as 
boxes. The images are in units of 10~^^ erg s~^ cm~^ . One pixel is ^150pc across. 
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Figure 19. Map of the Ha relative velocity for NGC 6166 
in Abell 2199. The increasing and decreasing velocities in the 
Eastern part of the emission may be explained if the low velocity 
blobs are closer and therefore subject to less extinction. The Ha 
flux is overlain as contours. The scale is in units of km s~^ . One 
pixel is ^150pc across. 



effect since no abs orption is seen in the integrated spectra 
from the outskirts (|Edwardsll2007[ ). 

The only region for which an H/3 emission strength is 
measurable is Region 8, which shows an Ha/H/3 ratio of 
2.4±0.4. This is within the theoretical lower limit of 2.8 to 
within our measureme nt errors. The ratio is m uch lower than 
the integrated value of lCrawford et al.l l| 19991 ). Ha/H/3 = 9.3. 
However, as Ha is strong in several regions and H/3 is only 
present in Region 8, the integrated ratio is naturally much 
larger. Nevertheless H/3 is clearly detected in Region 8 which 
suggests that the extinction level is lower in this region. It is 
therefore reasonable to suggest that Region 2 may be closer 
to the dusty torus of an AGN. It may also be that Region 8 
has a population of young stars and that H/3 absorption 
could be lowering this ratio. 



3.9 Cygnus-A 

Morphology - Figure [501 presents the continuum sub- 
tracted Ha, H/3, [N ll] A 6584, [O ill] A 5007 emission, and 
the continuum near Ha and near H/3. In the line emission 
images, it is the NW and the central components (Regions 1 
and 2) that have the highest flux. In the image of the con- 
tinuum near Ha (between 6380 and 6430 A), the brightest 



emission is seen in four blobs: one central, one East of cen- 
ter, one to the SE and another to the NW. The image of the 
continuum near H/3 (between 4720 and 4830 A) shows bright 
emission in the SE and NW blobs, but the central peak is 
not s een, likely due to intense dust extinction l|jackson et al.l 
I1998D . 

Kinematics - The relative velocity maps of the Ha 
and [O III] A 5007 emission lines are shown in Figure 1211 
The relative velocity map of the [N ll] A 6584 emission is 
not shown as it harbours very similar morphology and mag- 
nitude as to that for the Ha relative velocity map. Large 
scale velocity gradients of ±200kms~^ are seen for Ha, and 
of ±100 kms"^ for [O ill] A 5007, both in the SW and NE 
direction. On this scale, the velocity gradients in both lines 
are reminiscent of rotation. 

The FWHM of the Ha emission line is also shown in Fig- 
ure [21] It reveals large widths corresponding to the central 
emission peak, the loca tion of the central radio point source 
l|Tadhunter et al.ll2003D . There is also an area of large width 
that extends west of this peak (that is, west of Region 2). 
This second area extends to Region 8, the direction of the 
previously detected jet. The large width could be caused by 
disturbance from the jet. 

The morphology and velocity s een in these maps ar e 
consistent with the picture given by Ijackson et all (l998l ). 
Figure [22] summarizes their geometry and highlights an 
opening cone that passes from the NW side (Regions 1 and 
2) of the image through the central blob to the SW side 
of the image. These authors also point out blue condensa- 
tions (Region 6) and the jet that emanates from the center 
westwards, passing along south side of the NW blob. Our 
largest velocity gradients are located perpendicular to the 
op ening cone as well as to the direction of the jet identified 
by I Jackson et all l|l998D . 

Emission Diagnostics - Figure [3] shows the red and 
blue spectra for the central Region 1. As suspected, strong 
[O III] A 5007 and [N ll] A 6584 relative to H/3 and Ha. Fig- 
ure [4] shows that most pixels have [N ll] A 6584/Ha ratios 
higher than ^ 1.4, placing them on the AGN side of the 
BPT diagram. The highest ratios are in two clumps almost 
directly South and North of the center of the image, however 
there appears to be no relationship between the Ha lumi- 
nosity and the value of the line ratio. The lowest line ratios 
are East of the center, at Region 6. This is the location of 
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Figure 20. Cygnus-A. Top i?oui; From left to right, the continuum subtracted Ha emission flux, the continuum subtracted HfS emission 
flux, and the continuum subtracted [Nil] A 6584 emission flux. The red emission lines show very similar morphology, however, the [Nil] 
emission is much stronger, we are looking at gas clearly ionized by an AGN. Bottom Row: From left to right, the continuum subtracted 
[O in] X 5007 emission flux, the continuum near the Ha emission line and the continuum near the Hf) emission line. The central peak 
seen in the red continuum is not apparent in the blue continuum, as it is probably hidden by dust extinction. The Ha emission is overlain 



as contours. The regions used in the analysis are represented as boxes. The image is in units of 10 
^220 pc across. 
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Figure 21. Map of the line kinematics in Cygnus-A. From left to right - the Ha relative velocity, the [O ill] A 5007 relative 
velocity, and the Ha FWHM. Large scale velocity gradients are due to rotation. The largest values of FWHM are found at the center 
(Region 2) near the central black hole and also between Regions 3 and 8, near the location of the jet. The Ha emission overlain as 
contours. Notice that the [O ill] image needs to be shifted North a few pixels to align properly to the Ha image, and in the central panel 
the white pixels at the top represent regions of no data. The scale is in units of kms~^ . One pixel is ^220pc across. 
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Figure 22. Cygnus-A Ge o metr y. This image is adapted from 
Figure 4 ™ \jackson et al\ tl 99 j ) and shows the geometry of 
Cygnus-A. The position and direction of the opening cone and 
jet are shown on the a figure of the [ OI] line emission. The num- 
bers refer to the regions studied in this paper. The scale is in 
arbitrary flux units. 

bright, blue, condensed clumps of continuum emission (ob- 
se rved using the F622W HST filter) previously discussed 
in Ijackson et al.l ()l998l l . The spatial resolution here is not 
as high as for HST, and neighbouring pixels may be wash- 
ing out signs of young stars. Indeed, Ijackson et alj (|l998l ) 
are unable to identify any emission of [O ill] A 5007 or H/3 
within these blue condensations. 

This is the only BCG for which we can constrain the 
emission mechanism using all four diagnostic lines. The 
[O III] A 5007/H/3 map (not shown) is affected by the poorer 
quality of the H/3 image, but shows no clear structure im- 
mediately apparent in the map. However, all of the central 
pixels show ratios of [O ill] A 5007/H/3 >8, disclosing the na- 
ture of this AGN to be a Seyfert. The BPT diagram shown in 
Figure [53] classifies virtually every pixel as ionized by AGN 
emission, and there is a smooth transition from pixels at the 
Seyfert corner down to the LINER side. This agrees with the 
well established classification of Cygnus-A being a narrow 
line radio galaxy. Two additional diagnostics are available 
for Cygnus-A, the ratio of [O l] A 6300/[O ill] A 5007 and 
[O III] A 5007/HQf. The former is a good tracer of the ion- 
ization parameter (q, as defined in Kewley & Dopita 2002; 
q=uxc, a measure of the number of ionizing photons per 
atom at the boundary layer). Both lower and higher ion- 
ization states are present. The [O ill] A 5007/Hq: ratio also 
shows a region of higher ionization along the edge of emis- 
sion, North of the NW blob (Region 1), which was noted in 
the observations of ljackson et al.l l|l998h . 



4 DISCUSSION 

4.1 Overview of Results 

We have investigated the emission line morphologies, dy- 
namics and ionization state in a sample of nine BCGs. 
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Figure 23. BPT diagnostic diagram for Cygnus-A. Indi- 
vidual pixels are plotted as open circles and pixels within the ana- 
lyzed regions are plotted as filled circles. Almost every pixel shows 
exhibits the well known Seyfert status of Cygnus-A. This figure 
uses the definitions of H ll region, Composite and AGN ionizing 
sources from Kauffmann et a l. (2003). The desi gnations between 
LINER and Seyfert are from lOsterbrock & Ferlano, C2006,) . 

Condensed (Abell 2052, Cygnus-A), filamentary and patchy 
(AbeU 1060, AbeU 1668, and MKW3s), and extended 
(Abell 1204, Abell 2199) morphologies are observed and the 
Ha and [N ll] emitting gas usually follows the same morphol- 
ogy. Also of note are two cases, Abell 1060 (a non-cooling 
fiow) and Ophiuchus (at best a low level cooling flow), where 
the line emitting gas morphology echos that of prominent 
dust features seen in the optical images. By using BPT dia- 
grams to diagnose the ionization mechanism, we found that 
in two out of seven emitting BCGs, hot stars are likely 
perpetrators - NGC 3311 in the non-cooling fiow cluster 
AbeU 1060 and the BCG in cooling fiow cluster AbeU 1204. 
In these cases, SFRs and ages for the young populations 
are derived. However, the presence of an optical AGN, in 
both cooling flow as well as non-cooling fiow systems, was 
strong. Two of the nine BCGs observed had no emission 
lines. This includes the BCG of AbeU 1651 and the BCG 
Ophiuchus. Unexpectedly in Ophiuchus, it is another source 
within the cluster (Object B), rather than the BCG that 
shows line emission. In most cases, the relative velocities 
are in the neighborhood of ±100 - 200kms~^, though some 
higher relative velocities indicative of large infiows (Object B 
in Ophiuchus) or outfiows (MKW3s) are observed. 

The properties and main results for each cluster BCG 
are highlighted in Tableland discussed briefiy below. 

• A bell 1060: Previou s detections of star formation in this 
dusty l|Laine et al.ll2003l ') non-cooli ng fiow (iHavakawa et al.l 
120061 ) BCG have been described bv lVasterberg et all l|l99ll ). 
We characterize several regions of star formation, note that 
the morphology of the star forming region follows that of the 
dust, and present a smooth velocity gradient reminiscent of 
rotation. 

• Abell 1204: This cooling fiow cluster l|Bauer et al.ll2005l ) 
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has the highest redshift of our sample (z — 0.1706) and 
thus a larger extent of the cD galaxy is viewed. Regions of 
ionization due to AGN are separated from those indicative 
of a young stellar population. The latter are found further 
from the center in a plume that extends towards a chain of 
smaller cluster galaxies. The relative velocities are modest 
(±200 kms~^) and show no signs of rotation. An interaction 
with the nearby galaxies may be important. 

• Abell 1668: This is a non-cooling flow cluster 
llSalome fc Combed l2003l l. yet strong lines are present. 
These appear in a filamentary distribution with line ratios 
signifying ionization from an AGN. The relative velocity of 
the line emitting gas shows that it is not at rest with the 
underlying cD. 

• Ophiuchus: This cluster was recently observed to have a 
cool core, even though the X- ray tem perature of the cluster 
itself is very hot at 9-10 keV (iFuiita e t al. 2008). The BCG 
shows no emission lines. However, emission lines are seen in 
Object B, at a projected distance of ~2kpc away from the 
center of the BCG, which is also at the position of the X- 
ray centroid. Object B is also cospatial with a dust feature 
seen on the acquisition image, and the emission lines have 
a relative velocity of -1-750 km s~^ with respect to the BCG. 
We attribute this to infall onto the BCG. 

• MKW3s: This cooling flow cluster (Kaastra et al. 2004) 
has previously been observ ed to have a UV excess and at - 
tribute d star formation by McNamara fc O'Connelll (Il989l ) 
and by iHicks fc Mushotzkvl (|2005i r We find the emission 
lines in the GMOS IFU image to be filamentary in morphol- 
ogy, yet well described by ionization from an AGN. The lines 
are blueshifted by -1-560 km s~^ with respect to the BCG, 
suggesting an outflow. 

• Abell 1651; This cooling flow cluster l|Whitell2000l ) BCG 
shows no emission lines. 

• Abell 2052: This coolin g flow cluster (jBlanton et al.l 
I2OO3I : iKaastra fc et al] l2004l) has patc hy dust in the cen- 
ter of the BCG (iLaine et al.ll2003h and lHicks fc Mushotzkvl 
(|2005l ) and iBlanton et all l|20oi r have deduced star forma- 
tion in the BCG from excess UV-IR and U-band continuum 
emission. We find that the Seyfert signature overwhelms any 
star formation in the central few arcseconds of the BCG. The 
morphology of the emission lines here is point-like, in that it 
shows no asymmetric or extended features. The relative ve- 
locities vary on the scale of ±250 km s~^ , but are not smooth 
enough to warrant a classiflcation of rotation^ 

• A bell 2199: Another cooling flow jjohnstone et al.] 
|2002| ) BCG with pre vious calculations of star form ation 
based on UV excess (jMcNamara fc O'ConneUl Il989^ . We 
find the AGN signature in the emission lines is dominant 
throughout the extent of the OASIS image. The lines are 
brightest in the center, but exhibit an extended morphology 
towards the East. The relative velocities are similar to those 
seen in Abell 2199, but again the gradient is not regular 
enough to signify rotation. 

• Cygnus-A: T his poor cooling flow cluster 



Reynolds fc Fabiad Il996l) 



a weU studied AGN 



Jackson et al.lll998l : iTadhunter et al ]|l994l. l2003h. We find 



a mor phology consistent with the images of Jackson et al.l 
l|l998h . line ratios telling of an AGN, and a velocity gradient 
that suggests rotation. The FWHM of the emission lines 
traces the direction of the jet, as it becomes wider for the 
lines in that region. 



Certainly gas exists in BCGs and is excited, and the 
observations above show that many mechanisms are at play. 
For instance, signs of emission by hot stars are present in 
both cooling fiow and non-cooling fiow systems, and AGN- 
ionized gas is also present in both cooling flow and non- 
cooling flow systems. 

4.2 Scenario 



IWilman et al.l l|2006h . who conducted a similar study, find 
little variation of the emission line ratios across the emission 
nebulae in their high redshift and very luminous sample of 
4 cooling flow BCGs, implying a uniform ionization state. 
They suggest the following single scenario to explain the 
line emission in all their galaxies: an interaction of smaller 
cluster galaxies triggers starbursts in cold gas reservoirs, pre- 
sumably deposited from the cooling flow. This is supported 
by their observation that different ionization states of the 
gas vary little spatially, implying a single ionization source 
for the Ha emission. Intriguingly, the current observations 
for the modest cooling flow case of Abell 1204 support this 
idea (this is the only cluster in this sample where z > 0.1). 
In this case, the Ha, [N 11] and [S 11] emitting gas all share 
a similar morpho l ogy. B ut, in contradiction to the findings 
of IWilman et al.l (|2006l ). we do find a difference in line ra- 
tions with respect to position and Ha strength with the 
lower line luminosity pixels showing more star- forming like 
activity (consistent with results of Hatch et al. 2007). We 
calculate a total SFR of ~7M0 yr~ ^ for this BCG, wh ich is 
close to the infrared derived SFR of lO'Dea et~aLl (120081 ) and 
not inconsistent with the mass deposition rate as not all the 
molecular gas will convert to stars and as the calculations of 
MDRs based on the absence of the coolest gas are an order 
of magnitude below the rat es based on the classical cooling 
flow scenario - 50 yr""*^ (lO'Dea et al..,2008) co mpared to 
675 Mq yr" ' 



- ou ivip ) yr (| 
l|White|[2OO0l ). 



IWilman et al 



(|2006l ) also rea- 



soned that the Ha and CO gas are related as they share 
the same kinematics, and it is this CO gas, which is sub- 
sequently disturbed by the passing of nearby neighbors or 
by an AGN, which will emit the Ha line. The observations 
presented here for Abell 1204 show no hard evidence for an 
interaction with the nearby neighbors, but, it is tempting 
to speculate on the likelihood of an interaction with nearby 
galaxies seen on the acquisition image. It would be interest- 
ing to search for molecular gas in this BCG. This is the only 
cluster i n the sample present ed in this paper which agrees 
with the lWilman et al.l (120061 ) hypothesis. 

However, the overall "mixed bag" of ionization sce- 
narios seen throughout most of this sample i s even more 
appea ling when put into the context of the iHatch et al.l 
(120071) observation s. In o pp osition to th e scena rio put forth 
by IWilman et al] (|2006| ). iHatch et all (|2007l ) find cases 
where the line emission properties in cooling flow BCGs 
suggest motions from strong AGN or starburst driven 
outflows (Abell 2390, Abell 1068), from galaxy passbys 
(2 A 0335±096), and from rotation (Abell 262). They find 
that the ionization state is not uniform, and do not con- 
clude that one scenario can account for the emission lines 
seen in cooling flow BCGs, the results presented here are 
consistent with such conclusions. 

One caveat of our study is the small sample size, 
nonetheless, no difference in the ionization mechanism has 
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been found that clearly separates the cooling flow and non- 
cooling flow cluster BCGs. Neither has a consistent picture 
been developed to explain the origin of the line emission 
throughout the sample. Although, the emission line charac- 
teristics are consistent within each BCG and most systems 
show a hard ionizing source prominent throughout the cen- 
tral few arcseconds, usually well described by LINER emis- 
sion line ratios. 



As found by seve r al other authors ijWilman et al. I I2OO6I : 
iBohringer et al.ll2002l : iDonahue et al.ll2000l ). this data does 
not support a simple picture in which X-ray gas cools into 
molecular clouds subsequently forming stars. This is empha- 
sized by the variation in morphology of the ionized nebula in 
these systems. The disturbed morphology of the non-AGN 
ionized gas of Abell 1204 extends towards the direction of 
several smaller galaxies, suggestive of an interaction with 
nearby companions. The spectrum of MKW3s is also in- 
teresting, with the Ha emission line shifted with respected 
to the underlying spectrum. Although for the majority of 
cases the line emission is stronger in the central regions of 
the BCG itself, in Ophiuchus at least, the line emission is 
localized North of the BCG center, at a distance of '~2kpc. 
iDonahue et al.l (|2007l ) find the emission lines in the cooling 
flow cluster 2A0335+096 are associated with the BCG, as 
well as a companion. In the case of 2A0335-I-096, the emis- 
sion lines are described as dusty, and this is also the case in 
Ophiuchus and in NGC 3311 of Abell 1060 where the fine 
emission is constrained by what appears on the acquisition 
images to be strong dust features. 

There are of course many possible outcomes for any 
molecular gas that drops out of the cooling flow onto the 
BCG. In general, the observations from the dataset pre- 
sented here, of an overwhelming influence of AGN signatures 
in most of the sample support the scenario currently put 
forth for the nature of cooling gas in X-ray clusters where 
AGN feedback is important. That is, part of the material 
condensing out of the cooling X-ray gas finds itself in the 
form of molecular reservoirs at the centers of the BCG, and 
part rains onto the central black hole. This could trigger 
an outburst from the AGN which is then hypothesized to 
reheat the cluster enough to prevent any further cooling. 
The AGN could simultaneously ionize the molecular gas de- 
posits and infiuence starbursts. The molecular gas deposits 
would also be subject to fiybys from companions, another 
mechanism which could trigger a starburst. This complex 
scenario does explain the observations of line emission re- 
sulting from the ionization of hot stars and AGN activity 
in cooling flows. It could also explain observations of AGN- 
ionized line emission in non-cooling flow, as the line emission 
could be triggered after the cooling cluster gas has been re- 
heated by the AGN. The non-cooling flow BCG NGC 3311 
of Abell 1060 ma y also be a part of the nominal ~15% of line 
emitting BCGs I Edwards et al.lf2007l ) as it is a non-cooling 
flow and the emission lines show no evidence for AGN ion- 
ization, and thus perhaps the star formation is due to a 
completely different process. It is however, somewhat more 
difficult to explain the observations of Ophiuchus, where the 
line emission is not localized in the BCG. In a future paper, 
we will combine X-ray and radio observations with our non- 
detection of line emission from Ophiuchus and Abell 1651 to 
make a census of the heating and cooling in those systems. 



5 CONCLUSIONS 

The observations from integral field spectroscopy have re- 
vealed the complex nature of the line emission in these galax- 
ies, for which star formation, AGN, or both are important 
in several systems, but for which no consistency of the emis- 
sion line characteristics and the cluster properties is seen. 
The line emission is point-like, filamentary, or extended; it 
exists in regions that are plagued with large patches of dust, 
and those that are relatively dust free. The AGN ionizes all 
of the gas within the image in some cases, whereas regions of 
star formation are discernible in others. Outfiows are seen, as 
well as gas that has bulk rotations. And the emission could 
be associated with galaxy interactions in some, but not all 
cases. Each case with emission lines in a cooling fiow shows 
regions where the emission lines have AGN-like ratios. 

The process that is fueling the nominal emission in the 
non-cooling flows and the other bright galaxies at the cen- 
ter of a cluster is still unknown. However, in cooling fiows, 
the existence of an AGN is correlated with the presence 
of emission lines. As seen for example i n Abell 1204 (this 
work), 2A0335-K096 (|Hatch et al.l l200it ). and AbeU 2052 
(|Blanton et al.ll2003l '). the ionization mechanism can change 
with the radial distance from the center, showing lines ion- 
ized by a hard source at the center, and those ionized by hot 
stars further out. 

This work supports the current ideas on the impor- 
tant role of a massive black hole in these galaxies, as much 
of the optical line emission is dominated by gas ionized 
by an AGN. The black hole may a ct as a sink for the 
coolin g gas ([Pizzolato fc " Sokedl2005h . may trigger a star- 
burst l|Wilman et al.l 20061 ). and its energy output may pre- 



, a y 

vent further cooling of the cluster g as (Best et ahl lOOSal : 
ICroton fc et al.ll200d : iMcNamara fc Nuls en 2007i). This last 
role is an important factor in our understanding of the build- 
ing of large galaxies. The observed galaxy luminosity func- 
tion undergoes a sharp cut-off at the high end. This cut-off 
can be explained with the inclusion of X-ray gas and AGN 
feedback, where the amount of cooling (and hence poten- 
tial star formation) is regulated by heating from the AGN 
(|Croton fc et al.ll2006l) . 
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